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DETECTING COLD GAS AT Z = 3 WITH ALMA AND SKA 
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ABSTRACT 

We forecast the abilities of the Atacama Large Millimeter/submillimeter Array (ALMA) and the Square Kilo- 
meter Array (SKA) to detect CO and HI emission lines in galaxies at redshift z = 3. A particular focus is set on 
MiUcy Way (MW) progenitors at z = 3 for their detection within 24 h constitutes a key science goal of ALMA. 
The analysis relies on a semi-analytic model, which permits the construction of a MW progenitor sample 
by backtracking the cosmic history of all simulated present-day galaxies similar to the real MW. Results: (i) 
ALMA can best observe a MW at z = 3 by looking at CO(3-2) emission. The probability of detecting a 
random model MW at 3-cr in 24 h using 75 km s"' channels is roughly 50%, and these odds can be increased 
by co-adding the CO(3-2) and CO(4-3) lines. These lines fall into ALMA band 3, which therefore represents 
the optimal choice towards MW detections at z - 3. (ii) Higher CO transitions contained in the ALMA 
bands > 6 will be invisible, unless the considered MW progenitor coincidentally hosts a major starburst or an 
active black hole, (iii) The high-frequency array of SKA, fitted with 28.8 GHz receivers, would be a powerful 
instrument for observing CO(l-O) at z = 3, able to detect nearly all simulated MWs in 24 h. (iv) HI detections 
in MWs at z = 3 using the low-frequency array of SKA will be impossible in any reasonable observing time, 
(v) SKA will nonetheless be a supreme HI survey instrument through its enormous instantaneous field-of-view 
(FoV). A one year pointed HI survey with an assumed FoV of 410 deg^ would reveal at least 10^ galaxies 
at z = 2.95-3.05. (vi) If the positions and redshifts of those galaxies are known from an optical/infrared 
spectroscopic survey, stacking allows the detection of HI at z = 3 in less than 24 h. 

Subject headings: Galaxy; evolution - galaxies: ISM - radio lines: ISM - cosmology: theory 



1. INTRODUCTION 

Detecting cold gas in ordinary distant galaxies is a 
paramount challenge in modern astronomy. It will be ad- 
dressed by two revolutionary future telescopes, the Ata- 
cama Large Millimeter/submillimeter Array (ALMA) and the 
Square Kilometer Array (SKA), whose prospects have been 
stimulating much interest for cold gas at high redshift (z). 

Here, "cold gas" refers to all gas cold enough to be neu- 
tral. Such gas consists of molecules and atoms and it domi- 
nates the interstellar medium (ISM) in galaxies, mainly in the 
form of hydrogen and helium with a mass ratio close to 3:L 
The hydrogen is called HI when atomic and H2 when molec- 
ular Cold gas owes its astrophysical importance to its parti- 
cles moving slow enough for the formation of gravitationally 
self-bound structures on sub-galactic scales. Cold gas thus 
plays a primeval role in the formation of galaxies and stars 
(iKauffmann et al.ll2006t [Gunawardhana et al.ll201 ll) . Yet, the 
cosmic history of these processes remains puzzling, since ex- 
trapolations from available local observations to the past are 
complicated by the evolving conditions of the Universe, such 
as density, structure, and chemical composition. Hence, it is 
still unclear as to when the first galaxies formed, whether their 
cold gas was mostly atomic or molecular, how the first stars 
were born, and when the ISM was sufficiently chemically en- 
riched to form planets - to name but a few issues calling for 
detailed observations of distant cold gas. 

The difficulty of such observations arises from the extreme 
faintness of the cold gas tracers. Unless exposed to exciting 
radiation, cold gas only glimmers in narrow emission lines 
at infrared, (sub)millimeter, and radio frequencies. Molec- 
ular gas is most typically found via the J ■ 115 GHz (rest- 



frame) lines of the rotational transition s J ~ ^ of the 

'^C'^O molecule (hereafter CO, e. g. iTacconi et all 120 1 Ot 
lDaddietaLll2010l) . Atomic gas is detecte d via the "HI hne" 
or "21 c m Hne" at L420GHz rest-frame (IZwaan et al.ll2005t 
iMartin et al. 2010). To appreciate the difficulty of cold gas de- 
tections, note that the bolometric luminosity of a single bright 
star, such as Rigel (J3 Orionis), is 1000-times higher than the 
HI line power of the entire Milky Way (MW). Mainly for this 
reason, no HI emission has yet been seen at z > 1, while stel- 
lar light detections at such re dshifts are now an ob servational 
standard (e.g. Szomoru et al. 2 01 UlSu et al.ll201 ll) . 

The unprecedented sensitivity of ALMA and SKA in the 
(sub)millimeter and radio spectrum will ease the detection of 
CO and HI emission lines at high z. In fact, SKA was origi- 
nally conceived as a pure Hl-telescope (Wilkinson 1991), and 
the first science goal of ALMA is to "detect spectral line emis- 
sion from CO or CII in a normal galaxy like the MW at a 
redshift of z = 3, in less than 24 h of observation"^. In prepa- 
ration for ALMA and SKA, predictions of their findings are 
needed to optimize the telescope designs, to outline initial sur- 
vey strategies, and to ensure an unbiased check of our current 
theories against future observations. 

This paper illustratively predicts the abilities of ALMA and 
SKA to detect CO and HI emission lines fr om galaxies at z - 
3, base d on a semi-analytic galaxy model (lObreschkow et al] 
l2009ah . Of particular interest is the detection of lines in 
"galaxies like the MW at z = 3" - the first ALMA science 
goal. It is necessary to clarify whether this means a galaxy 
identical to the MW placed at a cosmological distance cor- 
responding to z = 3, or rather a plausible MW progenitor at 

' all science goals at http://almascience.eso.org/about-alma/full-alma 
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z - 1- The former interpretation is more common, but we 
here adopt the latter for it is perhaps more sensible in a study 
ultimately dedicated to the understanding of our own origins. 
Yet, this interpretation complicates the predictions as they re- 
quire a model of the MW at a cosmic time corresponding to 
z = 3, i. e. 11 billion years back in time. 

Section|2]reviews the S^-SAX simulatiorQ. Section [3] stud- 
ies the cosmic evolution of MW-type galaxies in S^-SAX with 
a focus on the CO and HI lines at z = 3. Section |4] sum- 
marizes current specifications of ALMA and SKA, based on 
which Sections |5] and |6] predict the detectability of emission 
lines from MW-type galaxies at z = 3 and arbitrary galaxies 
at z = 3, respectively. Section|7]concludes the paper. 

2. THE S^^-SAX SIMULATION 

Our analysis relies on S^'-SAX (lObreschkow et al.ll2009ah . 
a computer model of neutral atomic (HI) and molecular (H2) 
hydrogen i n galaxies^, whic h builds on the Millennium simu- 
lation ( Spr ingel et al.l l2005l) . The Millennium simulation is 
a gravitational A^-body simulation of about 10'° dark mat- 
ter particles in a cubic comoving volume of (SOO/ZiMpc)-'. 
It models the formation of cosmic structure down to galaxy 
haloes as low in mass as those of the Small Magellanic Cloud 
(SMC), while tracking features as large as the Baryon Acous- 
tic Oscillations (BAOs). The cosmological parameters of the 
Millennium simulation are h - 0.73, where the Hubble con- 
stant Hq = 100/1 km S^'MpC"', Qmatter - 0.25, t^bai-yon - 

0.045, Qa = 0.75, erg = 0.9. 

Through a post-p rocessing of the Millennium simulation, 
iDe Lucia &Blaizotl ([2007, see also ICroton et al.ir2006) stud- 
ied the evolution of idealized model-galaxies placed at the 
centers of the dark matter haloes. The global galaxy prop- 
erties, such as stellar mass, cold gas mass, and morphology, 
were evolved according to discrete, simplistic rules. This 
"semi-analytic" processing resulted in a catalog of evolving 
and merging galaxies. The number of galaxies at a cosmolog- 
ical time of 13.7 ■ 10^ yrs (i. e. today) is about 3 ■ 10^, and each 
of these galaxies has a well-defined history of growing and 
discretely merging progenitor galaxies that have been stored 
in 64 discrete cos mic time steps. 

lObreschkow et al.l (I2009al) applied an additional post- 
processing to the galaxies in the semi-analytic simulation by 
iDe Lucia & Blaizol (l2007h in order to subdivide their cold gas 
masses into HI, H2, and Helium. They also assigned real- 
istic radial distributions and v elocity profiles to the HI and 
H2 components. Subsequently. Obresch kow et al.l(l2009 b^ in- 
troduced a model to assign approximate CO line luminosi- 
ties to the molecular gas of each galaxy. This model relies 
on a single gas phase in thermal equilibrium with frequency- 
dependent optical depths, and it approximately accounts for 
the following mechanisms: (i) molecular gas is heated by 
starbursts, AGNs, and the redshift-dependent CMB; (ii) over- 
lapping clouds in dense and inclined galaxies cause CO 
self-shielding; (iii) in compact galaxy cores molecular gas 
transits from a clumpy to a smooth distribution; (iv) CO- 
luminosities are metallicity dependent; (v) CO-luminosities 
are always measured relative to the redshift-dependent CMB. 
The integrated CO and HI line luminosities were further ex- 
panded into frequency-dependent profiles - typically double- 
horn profiles - by applying mass models and random galaxy- 
inclina tions (sine-di stribution). The semi-analytic galaxy 
model (iDe Lucia & B laizot 2007J) with our additional proper- 

^ online access at http://s-cubed.physics.ox.ac.uk 



ties for HI, H2, and CO is called "S -SAX-Box" as a reminder 
that the simulated evolving galaxies are contained within the 
cubic volume (box) of the Millennium simulation. 

Given S^-SAX-Box w e then constructed a virtual sky 
(lObreschkow et al.1 l2009c) by mapping the Cartesian coordi- 
nates (x, y, z) of the simulated galaxies ont o apparent posi - 
tions (RA, Dec, z), using the method of Blai zot et al.l (120051) . 
Alongside this mapping, the intrinsic CO and HI luminosities 
of each galaxy were transformed into observable integrated 
line-fluxes. The resulting virtual sky simulation is called "S^*- 
SAX-Sky" in conti-ast to S^'-SAX-Box. The maximal field-of- 
view (FoV) of S^-SAX-Sky depends on the selected maximal 
redshift z. At z = 3 the FoV is approximately 37.2 deg^, corre- 
sponding to the comoving surface area of (500//? Mpc)^ of the 
Millennium box, which is large enough to suppress significant 
effects of cosmic variance. 

In this paper we are using both S^ -SAX-Box and S^-SAX- 
Sky. S^-SAX-Box contains the pointers needed to back- 
track the cosmic evolution of MW-type galaxies to z = 3 
(Section [3]), while S^-SAX-Sky provides the apparent posi- 
tions and line fluxes required to study the detectability of 
the simulated galaxies (Sections |5] and |6l). Throughout the 
whole paper, we assume that the continuum emission can 
be perfectly subtracted, such that the line emission can be 
studied independently. For other assumptions, limitations, 
and uncertai nties of the S^-SAX simula tion, please refer to 
Section 6 in 'Obreschko w et al.l (l2009al) and Section 6.2 in 
lObreschkow et al. (2009b). 

3. EVOLUTION OF SIMULATED MW-TYPE GALAXIES 

We shall now investigate the cosmic evolution of the CO 
and HI line signatures of the galaxies "like the MW" in the 
S^^-SAX-Box simulation (see Section|2l). 

3.1. Definition of simulated MW-type galaxies 

By definition, we call a model galaxy at z = a "MW-type", 
if its morphological type, derived from the bulge-to-disk ratio 
(see eq. 18 in Obreschkow et al. 2009a), is Sb-Sc, and if it 
matches the stellar mass Ms, the HI mass Mhi, the H2 mass 
Mh,, the HI half-mass radius r^^\ and the H2 half-mass ra- 
dius r^^^ of the MW, given in Tab.[T] within a factor 1.3. This 
factor approximately corresponds to the empirical uncertain- 
ties. According to this definition, the S^-SAX-Box simulation 
contains 1928 MW-type galaxies at z = 0. A simulated galaxy 
at redshift z > is called a "MW-type" galaxy or a "MW pro- 
genitor", if, at its particular redshift, it is the most massive 
progenitor of a MW-type galaxy at z = 0. 

3.2. Evolution of HI and H2 in simulated MW-type galaxies 

S"'-SAX-Box consists of 64 discrete cosmic time steps 
(Section|2]i. The cosmic evolution of any galaxy through those 
time steps can be extracted using a system of galaxy identi- 
fiers and progenitor-pointers that was alr eady installed in the 
underly ing semi-analytic g alaxy model (ICroton et al.l l2006l 
see also lSpringel et alj|2005i) . We here used those pointers to 
follow the cosmi c his tory of the 1928 present-day MW-type 
galaxies (Section lTTI ). 

The cosmic evolution of the sample averages of the masses 
and radii of the simulated MW progenitors is displayed in 
Fig.[Tl and the specific values at z = 3 have been summarized 
in Tab. [1] We emphasize that the sample size of 1928 MW- 
type galaxies at z = decreases monotonically with z, since 
the different evolution scenarios of the MW-type galaxies start 



Quantity Obs. atz = Sim. at z = 3 Ref. 



Virial mass Mvi,- [Mg] 


1.3±0.3- 10'2 


7 ^+0.9 
^"'-0.7 ' 


10" 


(a) 


Stellar mass Ms [Mq] 


5.0^1 ■ 10'" 


5 3+3' ■ 


10' 


(b) 


HI mass Mm [Mg] 


8.0+2 . iq9 


/hi • 0.7 


+1.1 . 10' 

-0.5 


(c) 


H2 mass Mht [Mq] 


3.5+1 ■ 10' 


2.8+;i . 


10' 


(d) 


HI half-mass rad. [kpc] 


15- 






(c) 


H2 half-mass rad. rj^"'' [kpc] 


7+] 






(d) 



TABLE 1 

Observed properties of the MW versus properties of the simulated 

MW-PROGENITORS AT Z = 3. ThE INDICATED RANGES ARE l-CT UNCERTAINTIES FOR 
THE OBSERVED VALUES AND RMS SCATTERS AROUND THE SAMPLE AVERAGES FOR T HE 
SIMULATED VALUES. ThE HI FUDGE FACTOR /hI IS EXPLAINED IN SeCTIOn IT?) 

. The observed values have been drawn from the following references: (a) 

McMillan (2011), (b) Flynn et al. (2006), (c) analytic fits to £ Hi(r) in 
Kalberla & DedesI 120081) , (d) I.H2 (r) in Tab. 3 in jSanders efal] fMl) . 

at different initial redshifts, depending on the respective dark 
matter distribution. In other words, with increasing z, the 
average MW properties displayed in Fig. [T] are increasingly 
biased towards evolutionary scenarios, which started particu- 
larly early in the history of the Universe. At z = 3 this is not 
an issue, since in 90% (i. e. 1731) of all simulation scenarios 
the MW formed before z = 3. However, only 39% (i. e. 755) 
of all scenarios has the MW forming before z - 1, and only 
12% (i. e. 234) of the scenarios before z = 8. 

Detailed physical interpretations of the cosmic evolution 
displaye d in Fig. [T] can be found in lObreschkow & Rawlingsl 
(l2009bl evolution of t he HI and H2 masses) and 
lObreschkow & Rawlingsl (1200 9a', evolution of the galaxy 
sizes). In brief, the radius of individual disk galaxies grows in 
the simulation with cosmic time approximately as (1 + z)"', 
consistent with optic al /infrared high-redsh if t observations 
jBouwens et al.1 I2004t iTruiillo et al] I2006t iBuitrago et all 
1200a) . This size evolution is reflected in the evolution of 
the HI and H2 radii (see Fig. [Tj bottom panel), and it is 
responsible for an increase in the pressure of the interstellar 
medium (ISM) with z. By virtue of the rel ation between the 
ISM pressure and the H- ) /HI ratio (e. g. lElmegreenI [19931: 
iBlitz & Rosolowskvl2006flLerov et al.ll2008h . flie H9/HI mass 
ratio therefore increases with z. This results in a roughly 
constant H2 mass for the simulated MW galaxies in the 
redshift range z = - 3 (see Fig. [Tj top panel), while the 
HI mass varies by a fact or 10 in the same redshift range (but 
see discussion in Section [T4] l. 

The simulated MW progenitors at z = 3 exhibit an average 
gas mass fraction (Mm + Mh,)/(Ms -i- Mhi + Mh,) of 40%, 
respectively 55% when correcting the HI masses as described 
in Section [374] These values lie an order of mag nitude above 
those found in today's massive spiral galaxies (iLerov et"al] 
120051) . and they are in good agreement with the average gas 
mass fraction of 44%, recently measured in typical massive 
star-forming galaxies at z ~ 2.3 (iTacconi et al.ll2010l) . 

3.3. CO and HI lines of simulated MW progenitors at z — 3 

Fig. 12] (data available onlin^ displays the emission line 
characteristics of our sample of 1731 simulated MW progen- 
itors at z = 3. The upper panel shows the sample distribu- 
tions of the frequency-integrated line fluxes S , while the lower 
panel represents th e distributions of the peak fl ux densities Sp. 
Using eq. (All) in lObreschkow et akl (l2009bl) the frequency- 

3 http://s-cubed.physics.ox.ac.uk/downloads/mw-at-z3.xls 
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Fig. 1 . — Cosmic evolution of the average properties of the simulated MW 
progenitors. The average values (solid lines) are arithmetic sample averages. 
The shaded regions represent the 0.5-0" scatter around the average values in 
the Gaussian approximation, i. e. those shaded regions contain about 40% 
of all events. Note that the number of galaxies, from which the averages 
and scatters were calculated, decreases with redshift, because the dift'erent 
evol ution scenarios for the MW start at dift'erent cosmic times (see Section 
13.21 . Also note that this fi gure shows the HI masses without the correction 
introduced in Section [3!4l 

integrated fluxes S (here in units of Wm"^) can be converted 
into velocity-integrated fluxes (e. g. in units of Jykms"'). 

The first conclusion from Fig.|2]is that the sample of simu- 
lated MW progenitors covers a wide range of fluxes for each 
individual emission line. In fact, the root-mean-square (RMS) 
scatter of the line fluxes varies between 0.5 and 1 dex for the 
different lines. For the CO transitions up to CO(5^) the sam- 
ple distributions are roughly Gaussian in log-space, reflecting 
the underlying sample scatter in the H2 mass. However, the 
higher-order CO transitions are skewed towards the high-flux 
end in the distribution. For example, the highest CO(10-9) 
fluxes in the sample lie nearly five orders of magnitude above 
the sample median. This non-Gaussian excess of high fluxes 
for the higher-order CO transitions reflects the relatively rare 
cases where the molecular gas is heated by a massive starburst 
or an AGN. In fact, all simulated MW scenarios occasionally 
undergo starbursts and phases of intense black hole accretion, 
but at any given cosmic time, such as the time correspond- 
ing to z = 3, only a minor fraction (« 5%) of all galaxies in 
the sample is subjected to such an exceptional source of heat. 
We also note that the sample distribution of HI fluxes is the 
only distribution with a non-Gaussian excess in the low flux 
regime. This feature is again attributed to occasional black 
hole activity, which, in the semi-analytic setup (.Croton et al] 



4 




Integrated line flux 5* |Wm" | 




10-* 10"' 10-'' 10"= 
Peak flux density Sp [Jy| 

Fig. 2. — Normalized sample distribution of the frequency-integrated line 
fluxes S (top) and peak flux densities Sp (bottom) for the 1731 simulated 
MW galaxies at z = 3. Dots represent median values. Note that the HI fluxes 
shown here are drawn from S'-SAX without the correction introduced in 
Section [X4l These data are accessible online^. 

12006'), results in a suppression of the cooling flow. Due to the 
large scatter and non-Gaussianity of the flux distribution in the 
sample, "average line fluxes" can be ambiguous or at worst 
meaningless. For example the average integrated CO(10-9) 
flux lies three orders of magnitude above the most probable 
integrated CO(10-9) flux. For this reason, we shall restrict 
our considerations to median values, where necessary. Those 
values have been marked as dots in Fig. |2] and are listed in 
Tab. E] (columns 16, 17) in Section|5] 

How do the simulated CO Spectral Energy Distributions 
(CO-SEDs) of our model MWs at z = 3 compare to real 
data? Recent observations of CO-SEDs in massive disk-like 
galaxies at z ~ 1 .5 found that those systems yield CO-SEDs 
similar to that of the MW, in contrast to the highly excited 
CO-SEDs typically obser ved in high-z submillimeter galaxies 
dDannerbauer et all2009t) . Fig.[3]demonstrates that these new 
data are roughly in line with the simulated MW progenitors, 
50% of which yield CO-SEDs only slightly more excited than 
the inner MW disk. These low-excitation MW progenitors 
dominate the predictions for ALMA and SKA in Sections |5] 
and|6] However, the simulation also predicts the existence of a 
minority of highly excited CO-SEDs (light gray in Fig.O cor- 
responding to the simulated MWs that underwent a starburst 
and/or an AGN at z = 3. 

3.4. The missing HI mass problem 

As discussed earUer dObreschkow & Rawlingsll2009bl) . the 
S^-SAX model misses a significant fraction of HI at z > 
compared to inferences from damped Lyman alpha systems 



■ BzK-21000(z=1.522) 
O MW inner disk 




~r 

3 4 5 6 7 8 9 
Upper rotational quantum number J 



Fig. 3. — Comparison of simulated CO-SEDs in MWs at z = 3 (shading) 
with observed data points of "normal" star forming galaxies: the inner MW 
disk (circles, Fixsen et al. 1999) and the massive galaxy BzK-2100() at z = 
1.522 (squares. Aravena et al. 2010; Dannerbauer et al. 2009). All SEDs are 
expressed in velocity-integrated fluxes, normalized to CO(2-l). The shadings 
respectively contain 50% of the simulated SEDs (dark gray), 70% (dark and 
mid-tone gray), and 90% (all gray tones). 

(DLAs). At z = 3, the global space density of HI in 
S^-SAX Hes a factor 5 below the DLA data Rec ent ef- 
forts to understand this difference dLagos et al.ll201 ll) claim 
that it can be widely explained by the limited mass reso- 
lution of the Millennium simulation, which defines a lower 
mass limit of Myjr ^ 10^"M q for the semi-analytic model of 
iDe Lucia & Blaizo3 (l2007h. Based on a Monte-Carlo extrap- 
olation to smaller galaxies iLagos etal] ( 1201 11) find that most 
of the HI gas at redshift z = 3 resides in galaxies that cannot 
be resolved in the Millennium simulation. This explanation 
of the missing HI mass in S^-SAX is further supported by 
high-res olution smoothed p article simulations of galaxy for- 
mation (iPontzen et alJl20()8h . which suggest that most DLAs 
are associated with small halo masses of M^h- = 10'^ - 1O"M0. 
While very plausible, these results remain uncertain because 
the H2/HI ratio of small high-z galaxies is poorly understood, 
in particular because their geometry is likely to deviate sig- 
nificantly from flat disk s. The cosmic evolution of metallicity 
and velocity dispersion (iForster Schreiber et alJl20b6) adds to 
this uncertainty. 

To address the many systematic uncertainties regarding 
HI at z = 3, we shall here consider two models: the raw S^- 
SAX model, which seems to underestimate the space density 
^Hi(z = 3); and a heuristic correction of the S^-SAX model, 
where all simulated HI masses are multiplied by a fudge fac- 
tor /hi = 5, matching Qhi(z = 3) inferred from DLAs. Phys- 
ically, /hi can be interpreted a correction containing the non- 
resolved satellites, as well as an Hl-rich non-disk component 
around each galaxy. All the HI detection predictions in Sec- 
tions |5] and |6| are provided for both the raw S^-SAX model 
(/hi =1) and the corrected one (/hi = 5). 

4. SPECIFICATIONS OF ALMA AND SKA 

This section outlines the provisional specifications of 
ALMA and SKA, needed for the predictions in Sections |5] 
and|6l Over the following paragraphs we present the physical 
concepts and assumptions behind the fundamental telescope 
parameters listed in Tab.|2]and the derived emission line spe- 
cific parameters listed in Tab. [3] 
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ALMA 


SKAi-LF 


SKA2-LF 


SKA2-MF 


SKAi-HF 


SKA2-HF 


Receiver type 


SFD 


AAS 


AAS 


AAS 


SFD 


SFD 


Diameter of dislies/stations D [m] 


12 


180 


180 


56 


15 


15 


Number of dishes/stations A'units 


50 


50 


250 


250 


125* 


1250' 


Number of instantaneous beams Wbeams 


1 


480 


4800 


4800 


1 


1 


RMS surface error of dishes 6 [mm] 


0.01 








0.5 


0.5 


Geometry factor 6p 


1 


eq. (|2) 


eq. (2) 


eq. (2) 


1 


1 


Con'elator quantization efficiency Eq 


0.95 


0.95 


0.95 


0.95 


0.95 


0.95 


Array efficiency fx 


0.90 


0.90 


0.90 


0.90 


0.90 


0.90 


Antenna efficiency 


eq.S 


0.90 


0.90 


0.90 


eq.d 


eq.(3) 


Nyquist-sampling frequency vo [MHz] 




115 


115 


800 






RMS baseline of compact configuration i?RMS [km] 


0.08 


100 


100 


100 


0.5* 


0.5* 


Receiver temperature Tree [K] 


Tab.|3] col. (|6) 


150 


150 


50 


30 


30 


Sky temperature Tsky [K] 


Tab. [3] col. (7) 


eq. (9) 


eq. 


eq. 


2.7 


2.7 



TABLE 2 

General specifications of ALMA and SKA. AAS stands for "aperture array station"; SFD for "single-feed dish". The values marked with ' correspond to 

THE CORE of SKA-HF. The full SKA-HF has twice the number of dishes and long baselines with Brms ~ 100 km. Core-only CO(l-O) observations are 
considered in this paper, since the high spatial resolution (i. e. low surface brightness sensitivity) of the full array would yield far less individual CO(l-O) 

detections at z = 3 (see Section|43). 



4.L Brief overview of ALMA 

ALMA is a reconfigurable array of 50 steerable single- 
feed dishes (SFDs). By exchanging the receivers, 10 differ- 
ent frequency bands, here called ALMA-1 to ALMA-10, can 
be reached. They collectively cover the whole atmospher- 
ically transparent parts of the spectrum between 31.3 GHz 
and 950 GHz. The window between 31.3 GHz and 84 GHz 
(ALMA-1 and ALMA-2) remains subject to future receiver 
development, and the availability of the window between 
163 GHz and 211 GHz (ALMA-5) is still uncertain, poten- 
tially impeding CO(6-5) and CO(7-6) observations at z = 3. 
The general ALMA specifications in Tab. |2] summarize the 
current online specification^ Before the completion of the 
array by 2013, the Large Millimeter Telescope (LMT) yields 
competitive sensitivities in the ALMA spectrum from ~ 80 to 
350 GHz. 

All redshifted CO emission lines at z = 3 considered in this 
paper are covered by the ALMA-bands except for the CO(l- 
0) line at 28.8 GHz, which is too low in frequency, and the 
CO(2-l) line at 57.6 GHz, which lies at the center of a ma- 
jor oxygen absorption band. For the remaining CO lines, the 
atmospheric transmissivity is remarkably high, such as illus- 
trated in Fig. nil for a low precipitable water vapour (PWV) of 
0.5 mm. In this paper, we adopt a slightly less optimistic value 
of PWV=1.0mm, which is realistic in the sense that lower 
(i. e. better) PWV-values have been measured over more than 
50% of the time over a year at the ALMA sitfl 

4.2. Brief overview of SKA 

SKA is only approximately specified and its concept might 
still change. Here we assume that SKA will be composed 
of three independent arrays with specificatio ns synthesiz- 
ing th ose de scribed by [P ewdnev et al. (201(|), iGarrett et al.l 
dSH), and ISchilizzi etal.l (l2007): the low-frequency ar- 
ray "SKA-LF", operating at 70MHz-450MHz (HI at z = 
19.3-2.2); the mid-frequency array "SKA-MF" at 400 MHz- 
1.4 GHz (HI at z = 2.6-0); and the high-frequency array 

^ http://science.nrao.edu/alma/specifications.shtml 

^ http://www.apex-telescope.org/sites/chajnantor/atmosphere 



"SKA-HF" at 1 GHz-30 GHz0 (e.g. CO( 1-0) line atz > 2.8). 
The latter is a fixed array of steerable SFDs, whereas SKA-LF 
and SKA-MF are fixed arrays of circular aperture array sta- 
tions (AASs) - a modern concept with no moving parts, cur- 
rently realized in the European Low-Frequency Array (LO- 
FAR). In this paper, we only consider SKA-LF (for the HI line 
at z = 3) and SKA-HF (for the CO(l-O) line at z = 3); but see 
the discussion in Section lT2l regarding potential uses of SKA- 
MF for HI detections at z ~ 3. 

SKA will be deployed in two phases referred to as SKAi 
and SKA2. The main differences between these phases are 
summarized in Tab. |2l In particular, SKA2-HF is assumed 
to have 10-times more SFDs than SKAi-HF, and SKA2- 
LF is assumed to have 5-times more AASs than SKAi-LF. 
Also note that in the current design, the mid-frequency ar- 
ray SKA-MF will only be added in phase-2. The comple- 
tion of SKA2 is not expected before 2022, i. e. around a 
decade after the completion of ALMA. Until then, several 
other telescopes and networ ks, summarized in Section 3 of 
iRawlings & Schilizzil ( 1201 Ih . will serve as technological and 
scientific SKA-pathfinders: the Australian SKA Pathfinder 
(ASKAP), the South African SKA Pathfinder (MeerKAT), 
the Westerbork Synthesis Radio Telescope (WSRT) upgraded 
with the phase array feed (APERTIF), the Murchison Wide- 
field AiTay (MWA), the Low-Frequency Array (LOFAR), the 
upgraded Multi-Element Radio Linked Interferometer Net- 
work (e-MERLIN), the electronic European VLBI Network 
(e-EVN), the European Pulsar Timing Array (EPTA), and the 
Five hundred meter Aperture Spherical Telescope (FAST). 
Also the following instruments will help preparing the way 
towards SKA: the Hydrogen Epoch of Reionization Array 
(HERA), the Extended Very Large Array (EVLA), the Giant 
Meterwave Radio Telescope (GMRT) with its new software 
correlator, and the Arecibo telescope via the Arecibo Legacy 
Fast ALFA Survey (ALFALFA). 



* The SKA design process currently uses 10 GHz as the required frequency 
upper hmit. However, it also tar gets 0.5 mm RM S surface accuracy to ensure 
high dynamic range imaging I G arrett et alJI2010D . Hence, the Ruze-equation 
(eq.[3) shows that the dishes should have good efficiencies up to 30 GHz. 
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Fig. 4. — CO emission line frequencies at z = 3 compared to the atmo- 
spheric transmissivity at the ALMA site (transmission function reprinted 
from https://almascience.nrao.edu/about-alma/aIma-site). All considered 
CO transitions except for CO(2-l) exhibit a high transmissivity, i. e. low 
atmospheric absorption. Note that the availability of bands 5 [CO(6-5) and 
CO(7-6)] is not yet confirmed. 



4.3. Point-source sensitivity 

The RMS noise cr (units proportional to Jy) of arrays of 
SFDs and AASs is approximated by (e. g. ,De Breuck. .2005) 



cr = 



2 Jc T^sys 



A fgeqfxea -y/npAvAf 



(1) 



where T^y^ = Tiec + T^ky is the system temperature that de- 
scribes the noise of the receiver (Tree) and the sky (Tsky) in 
the black body approximation; k is the Boltzmann constant; 
«p is the number of co-added polarizations (i. e. n^ - 2 if 
information on polarization is irrelevant), Av is the width of 
the frequency channels; At is the integrated observing time; 
and A is the physical surface area, that is literally the summed 
surface area of all SFDs or AASs. 

The dimensionless parameters eg, eq, e^, in eq. ([T]i are ef- 
ficiency terms, all between and 1, associated with the dom- 
inant causes of sensitivity loss. The dimensionless "geometry 
factor" Eg corrects A for geometrical projection effects and 
under-sampling effects in the case of AASs, while fg = 1 for 
dishes. Explicitly, 



mm 



cos 0, cos 0, 



for AASs, 



(2) 



where Q is the zenith-angle of the point-source, /lobs is the ob- 
serving wavelength, and /Iq = c/vq is the Nyquist-sampling 
wavelength of the array, eg = ■ fgi is the product of two 
interpretable factors: = cos corrects the collecting area 
for the linear projection of the wavefront onto the horizon- 
tal array; eg2 = min[l, /lobs//lo] min[l, /lob.s/C'^o cos 0)] ensures 
that the sensitivity is correctly reduced (cr is increased) in the 
case of a sparse (i. e. sub-Nyquist) sampling of the wavefront. 
The "correlator quantization efficiency" fq measures the noise 
level output by the correlator (software or hardware) com- 
pared to that of an ideal correlator; we here assume high-end 
correlators with eq - 0.95. The "array efficiency" ex mea- 
sures the losses due to the time needed to reconfigure the 
array and due the differential weighting app lied to the visi- 
bilities (i. e. tapering of the (u, v)-plane, see Holdawavl l 19981 ; 
l^n & Kogan 1999). Here we adopt the optimistic value of 
ex = 0.9 as we only consider fixed array configurations (the 
most compact configuration for ALMA and the permanent 
configuration for SKA), working in a low resolution mode 
close to the natural weighting of the visibilities (low tapering). 
Finally, the "antenna efficiency" (sometimes called "aperture 
efficiency") is the fraction of the electromagnetic energy 
transmitted from the collector to the receiv er For dish es, 
can be approximated by the Ruze-equation (lRuzelll952h 



ea = eo exp 



47t6 

^obs 



(3) 



where eo ~ 0.8 is the long wavelength maximum efficiency, 
S is the RMS value of the surface error distribution (values 
given in Tab. |2]l, and /lobs is the observed wavelength. 

We note that cr in eq. ([T]i is the RMS of the random compo- 
nent of the observing noise (in units proportional to Jy), mea- 
sured per frequency channel of width Av and per synthesized 
beam, i. e. per pixel of the synthesized sky image. Hence cr is 
the per-channel-noise of any source small enough to be fully 
contained within the synthesized beam. This feature applies 
in particular and always to point-sources. Therefore the sensi- 
tivity defined as the inverse of the noise cr is often referred to 
as "point-source sensitivity". For resolved sources, the sensi- 
tivity deteriorates, such as detailed in Section |43] 

4.4. Instantaneous field-of-view 

The full-width-half-maximum (FWHM) of the primary 
beam of an AAS or a SFD with diameter D is approximately 



FWHM K 1.22/lobs/£), 



(4) 



and therefore the field-of-view (FoV) of the primary beam is 

FoVbeam* ^^(/iobs/O)'. (5) 

AASs constantly collect electro-magnetic radiation from a 
large fraction of the hemisphere, typically covering about 
10 deg^ (half the hemisphere). The beams are "formed" 
through digital processing, and in principle the number of in- 
stantaneous beams A^beams is only limited by the processing 
power Given A^beams, the instantaneous FoV of AASs reads 



FoV = A^beams FoVbeam- 



(6) 



Note that we here assume that SKA-LF only uses digital beam 
forming. By contrast, the "pathfinder" instruments MWA and 
LOFAR also use analogue beam formers that cut down the 
FoV. 
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TABLE 3 

Considered emission lines of CO and HI with the corresponding telescope properties oe ALMA and SKA. Detailed descriptions of each column are provided 

IN Section ITtI 



4.5. Spatial resolution and sensitivity for resolved sources 

The spatial resolution y of telescope arrays corresponds to 
the beam width of the full array, which is also referred to as 
the "synthesized beam" to be distinguished from the "beam" 
of individual SFDs and AASs. We here approximate this res- 
olution by substituting D in eq. dU for the RMS length of the 
baselines Brms, 

7 = 1.22/lobs/BRMS, (7) 

although the precise resolution depends on the full baseline 
pattern, on the sky coordinates, and on the post-processing 
(i. e. weighting of the baselines). To check the resolution 
implied by eq. (|7]i for SKA-LF, we explicitly computed the 
point spread function (PSF) for a gen eric configuration of 
fifty AASs (as illustrated in Fig. 3 of Dewdnev et al.ll20Tol; 
station positions supplied by R. Millenaar, priv. comm.), cen- 
tered about (Ion 72°, lat -30°). The simulation consists of a 
complete track at 385 MHz, observing a source at a favorable 
declination of -60°, with flagging applied to scans below 30° 
elevation. The naturally-weighted restoring beam for such an 
observation corresponds to a resolution of 2.8" x 1.9" at a 
position angle of 95°, in close agreement with our derived an- 
gular resolution of 2.1" (Tab. [3] col.flSTl calculated via eq. (|7]i. 

In this paper, we focus on the detectability of individual 
galaxies. In order to obtain the highest sensitivity, it is there- 
fore desirable to chose the lowest possible spatial resolution. 
For ALMA, this is achieved by selecting the most compact 
array configuration (Brms - 0.08 km). This yields spatial res- 
olutions between 3.3" (CO(10-9) at z = 3) and 1 1" (CO(3-2) 
at z = 3), much larger than the average apparent H2 half- 
mass diameter of the MW disks at z = 3 of about 0.4", as 
derived from r^^^ in Tab. [T] using an angular diameter dis- 
tance of 1606 Mpc. No MW progenitor will hence be re- 
solved using the compact ALMA configuration. For SKA 
(Brms ~ 100 km) the situation is more subtle, since this in- 
strument is not reconfigurable. In the case of HI at z = 3, 
the resolution of 2.1" still exceeds the average HI half-mass 
diameter of the MW disks at z = 3 of about 1". No MW pro- 
genitor will hence be resolved. However, in the case of CO(l- 
0) at z - 3 imaged with SKA-HF, the resolution becomes 
as good as 0.03", such that a typical MW progenitor will 
be resolved in roughly 100 pixels (depending on inclination 



and evolution scenario). The associated order-of-magnitude 
loss in surface brightness sensitivity implies that many more 
sources will be picked up in core-only observations, where 
half of SKA's collecting area is sacrificed to the benefit of 
having all antennas within a core of 0.5 km radius. In this 
case the resolution drops to 5.2" and none of the MW at z = 3 
will be resolved in CO(l-O). In this paper we therefore as- 
sume that only the core of SKA-HF is used. 

Given these assumptions, none of the MW progenitors at 
z = 3 will be resolved in HI or CO emission. However, other 
galaxies z = 3, larger than the MW progenitors, may still be 
resolved. In this case, the sensitivity will be reduced relative 
to the point-source sensitivity of eq. ([T). In fact, if every pixel 
has an RMS noise level cr defined by eq. ([TJ, then a source 
extended over m > I pixels will be subjected to a Jy-noise 
equal to cr' = cr -\/m. Since real sources are not homogeneous, 
but rather exponentially fading disks, the definition of m is not 
obvious. Here we adopt the approximation 



where a is the angular half-mass diameter, measured along 
the major axis, of HI (for HI line) or H2 (for CO lines), and / 
is the galaxy inclination defined as the smaller angle between 
the line-of-sight and the galaxy's rotational axis. Our choice 
of using the the half-mass diameter rather than a larger diam- 
eter containing more of the gas mass relies on the assumption 
that the central concentration of gas can be exploited in clever 
algorithms for source extraction. 

4.6. Instantaneous bandwidth and spectral resolution 

ALMA's and SKA's limitations regarding the spectral reso- 
lution and the instantaneous bandwidth can be safely ignored 
within a study limited to the pure detection of extra-galactic 
emission lines in a narrow redshift range around z = 3 (Sec- 
tions |5] and |6j. For typical observing times (Af » 1/Av), 
the spectral resolution Av is mostly limited by the correlator 
performance. Both ALMA and SKA will be fitted with corre- 
lators allowing the selection of (frequency-dependent) spec- 
tral resolutions with an equivalent Doppler velocity far below 
1 km s- '. Th us the chosen velocity channels of 75 km s-' (see 
Section 14.71 ) never conflict with the spectral resolution limit. 
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As for the instantaneous spectral bandwidth (BW), ALMA 
(BW < 8 GHz) and SKA (BW < min[0.25 Vobs, 4 GHz]) are 
both able to cover a redshift range larger than Az - 0.1 at 
z = 3, which is the range that will be considered in Section|6l 
Even if the maximal instantaneous bandwidth were used, the 
implied highest spectral resolutions of ALMA (8192 chan- 
nels) and SKA (~ IC* channels) still provide channels much 
smaller than 75 km s"'. 

4.7. Perfonnance calculations for ALMA and SKA 

For all calculations hereafter frequency channels of an 
equivalent Doppler velocity of 75 km s ' are assumed. Since 
ordinary star forming galaxies at z « 2 show line widths of up 
to ~ 600k ms~' (FWHM wh en seen edge-on, Tacconi et al. 
120 lOl and iDaddi et al.l |2010j), channel widths larger than 
75 km s"^ might be beneficial for the pure detection of cold 
gas at z - 3. The predictions presented in this paper can be 
approximately rescaled to other channel widths w by multi- 
plying the signal-to-noise ratios n (see definition in Section|5]l 
by -\/vv'/(75kms"')- For example, a 3-cr detection (i. e. n = 3) 
with 75kms ' channels roughly corresponds to a 6-cr detec- 
tion (i. e. n - 6) with 300 km s ' channels. However, in prac- 
tise the signal-to-noise ratio of large channels would be low- 
ered by two mechanisms: the peaks in the line profile would 
get averaged-out, and a significant fraction of relatively face- 
on galaxies would have apparent line widths more narrow than 
the channel width. 

Tab. |3] displays the emission line specific performance of 
ALMA and SKA, as derived from the general telescope spec- 
ifications in Tab. |2] and the equations introduced in the pre- 
ceding part of this Section. The following list provides details 
and references for each column in Tab. |3] 

(1) Emission line identifier. 

(2) Telescope acronym. For SKA, the subscript speci- 
fies the construction phase, while LF (low frequency) 
and HE (high frequency) indicate the array type. For 
ALMA, the numbers indicate the frequency band. 

(3) Rest-frame frequency of the line center. 

(4) Observer-frame frequency of the line center, when ob- 
served at z = 3. 

(5) Channel width corresponding to an intrinsic velocity 
width (projected onto the line-of-sight) of 75kms 
calculated as Av - (75 km s^'/c) ■ Vobs- 

(6) Receiver temperature. For ALMA these values are 
drawn from the curre nt online specification fl ; for SKA 
they are adopted from lDewdnev et aTl (|2010|K 

(7) Sky temperature. For ALMA these values have been 
calculated as Tskv = Th^l + T.^tm(^ - ^) [see appendix 
eq. (1) in llshii et al.l l20Toll with a source tempera- 
ture T\, = Tqmb = 2.7 K (cosmic microwave back- 
ground), an atmosphere temperature ratm - 260 K, and 
atmospheric transparencies i] retrieved from ALMA's 
online'' 'Atmospheric transmission calculator" (for 
PVW - 1.0); for SKA-HF the sky temperature at 
28.8 GHz is considered equal to Tcmb; for SKA-LF 

' http://www.eso.org/sci/facilities/alma/system/frontend/ 



we follow (lBregmanll2000l see Fig. 1 in lDewdnev et al.l 
l20T0l) . 

r,ky=60K-(^bs/m)"5. (9) 

(8) System temperature Tjys - Tj^c + T'sky 

(9) "Geometry factor", representing the sampling effi- 
ciency of th e wa ve front in the case of aperture arrays 
(see Section l43] l. 

(10) "Correlator quantization efficiency", measuring the 
noise level added by the correlator (see Section l431 l. 

(11) "Array efficiency", representing the sensitivity losses 
due to tapering (see Section l431 l. 

(12) "Antenna efficiency" (also called "aperture efficiency"), 
representing the energy fraction actually transferred 
from the collector to the receiver (see Section l43] ). 

(13) Characteristic receiver noise cr given by eq. ([T]). 
Polarization information was ignored, i. e. «p - 2, and 
the observations are assumed close enough to the zenith 
that cos(ff) = 1 in eq. (O. 

(14) Instantaneous FoV calculated via eqs. (|5]|6]l. Note that 
the onUne ALMA specifications'* use FWHM a AohJD 
instead of eq. (HJi, giving a FoV 30% smaller than stated 
in Tab. [31 

(15) Spatial resolution calculated via eq. (|7|l. 

5. DETECTION OF MW-TYPE GALAXIES AT Z = 3 

Using the S^-SAX simulation (Section|2l) and the telescope 
properties of ALMA and SKA (Section HI we shall now in- 
vestigate the ability of these telescopes to detect the redshifted 
CO and HI emission lines of galaxies at z = 3. This section 
specifically addresses MW-type galaxies, while line detections 
in arbitrary galaxies will be considered in Section|6l 

From the wide range of possible observing goals at z = 3 
with ALMA and SKA we illustratively pick two questions: (i) 
What fraction of MW-type galaxies can be detected in each 
emission line at 3-cr and 10-cr significance in a 24 h pointed 
observation? (ii) What observing time is required to detect 
a random single MW-type galaxy with 50% chance? Here, 
"observing time" is defined as the integrated exposure time. 
It may, in practise, consist of multiple exposures spread over 
a period longer than the observing time itself. 

We remind that (Section O the 1731 simulated "MW-type 
galaxies at z = 3" are defined as the most massive progenitors 
at z = 3 (2.2 Gyrs after the Big Bang) of all simulated galaxies 
that have MW-like properties (see Tab.[TJ at z = (13.7Gyrs 
after the Big Bang). Therefore, the MW-type galaxies at z = 3 
represent a complete sample of MW progenitors within the 
semi-analytic model described in Section|2] 

By definition, a galaxy will be called "detected at n-cr" 
in a particular emission line, if the peak flux density Sp of 
that line lies n-times above the RMS noise cr of the obser- 
vation, i. e. the "signal-to-noise ratio" is n - Sp/cr. This is 
a conservative definition, since it makes no use of the spec- 
tral i nformation. Com bining multiple frequency channels 
(e. g. I Wang et aP 120061 for a related context) will undoubt- 
edly increase the significance of a detection, although such 
sophisticated techniques still require some development. 

T he fr action of MW progenitors (all non-resolved, see Sec- 
tion l4.5l l detected in 24 h at n-cr can be obtained by integrating 
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the normalized distribution of the peak flux densities Sp (lower 
panel of Fig. |2]l over Sp > n ■ cr, where cr - {cr yfKt) / V(24 h) 

with (cr VA?) drawn from Tab. [3] col. ( fT3] l. On the other hand, 
the observing time T required to detect 50% of the MW pro- 
genitors at n-cr is given by T = (cr VA?)^ rp- S^, where Sp is 
the median of the sample distribution of Sp. 

The results of these calculations are provided in Tab. ID The 
columns have been numbered such that Tab. |4] becomes an 
extension of Tab. [3] i. e. identical columns are given the same 
column number, while new columns are given a consecutive 
column number. These additional columns are: 

(16) Median value of log(5 / [W m"^] ) in the sample of 1 73 1 
simulated MW progenitors at z = 3, where S is the 
frequency-integrated line flux. 

(17) Median value of log(,s'p/Jy) in the sample of 1731 sim- 
ulated MW progenitors at z = 3, where Sp is the peak 
flux density of the emission line. 

(18) Fraction of simulated MW progenitors at z = 3 detected 
at 3-cr (or higher) in a 24 h observation. Values below 
0.1% are not resolved here, since they correspond to 
less than one simulated galaxy. 

(19) Same as col. ( fTST l but with lO-cr detection limit. 

(20) Exposure time required to detect a random simulated 
MW progenitor at z - 3 at 3-cr (or higher) with a chance 
of 50%. 

(21) Same as col. ( l20l i but with 10-cr detection limit. 

The results in Tab. |4] directly address the first ALMA sci- 
ence goal to detect spectral line emission from CO in a nor- 
mal galaxy like the MW at a redshift of z = 3, in less than 
24 h of observation. Col. ( fTSl l of Tab. |4] re veals that only the 
lines between CO(3-2) and CO(6-5) inclusive can make a 
serious contribution to the number of detections in the sense 
that each of those lines will be detected at a 3-cr level (or 
higher) in more than 10% of all model scenarios for the MW 
at z - 3. This fraction becomes maximal for the CO(3-2) 
line, contained in ALMA band 3. The same ALMA band 
also contains the redshifted CO(4-3) line, which has the sec- 
ond highest detection rate according to our predictions. We 
therefore conclude that ALMA band 3 can best respond to the 
first ALMA science goal. The odds for a MW detection with 
ALMA-3 can even be increased, assuming an ALMA correla- 
tor that allows the simultaneous observation of CO(3-2) and 
CO(4-3) at z = 3. In this case, the signal-to-noise ratio n of 
the co-added peak flux densities Sp ("single source stacking") 
can be approximated as 

'?p,CO(3-2) + ^p,CO(4-3) 

n = , (10) 

" CO(3-2) ^ " CO(4-3) 

We find that the odds of detecting a random MW progenitor 
at z = 3 (3-cr significance) are as high as 60%, if both the 
CO(3-2) and CO(4-3) line are used simultaneously. 

The higher order transitions including and above CO(7-6) 
will be virtually non-detectable. All those transitions require 
about or more than 1 yr (= 8.76 ■ lO-' h) of effective observing 
time to pick up a random MW progenitor with 50% chance 
at 3-cr significance. The few cases (< 10%) detected in 24 h, 
correspond to the objects in the pronounced wings of the flux 



distributions shown in Fig.|2] Those objects represent the sim- 
ulation scenarios where the MW underwent a massive star- 
burst or significant black hole activity exactly at a co smolog- 
ical time corresponding to z = 3 (see Section 13.3b . Those 
cases are therefore outside the scope of the first ALMA goal 
to detect a "normal" galaxy like the MW. 

CO(l-O) observations with the core of SKA-HF - should it 
reach 28.8 GHz - will be very powerful. SKA2-HF can detect 
more than 90% of all model MWs at z = 3 at 3-cr signifi- 
cance (or more) in 24 h. This high detection rate directly re- 
sults from SKA's huge collecting area, which widely exceeds 
that of ALMA, even when accounting for the relatively small 
aperture efficiency of - 0.56 at 28.8 GHz. 

On the other hand, SKA-LF will be virtually unable to de- 
tect HI emission from a MW progenitor at z = 3, as can be 
seen from the very long observing times (» 1 yr) required to 
pick up 50% of all MW progenitors in the simulation. Quite 
surprisingly, the detection rate of general galaxies in a blind 
1 yr HI survey with SKA2-LF will nonetheless be higher than 
that of any conc eivab le CO line survey with ALMA and SKA- 
HF (see Section l672] i. 

6. DETECTION OF ARBITRARY GALAXIES AT Z = 3 

This section expands the scope of Section |5] towards line 
detections in arbitra ry galax ies at z = 3 using ALMA and 
SKA. The Sections 16.1116.31 successively address three se- 
lected questions: (i) How many galaxies per unit sky area 
and redshift will be detected at 3-cr and 10-cr in a 24 h single 
pointing? (ii) How many galaxies in a redshift range Az = 0.1 
around z = 3 will be detected during a 1 yr survey? (iii) What 
is the significance of a stacked signal obtained in a 24 h ob- 
servation of a redshift range Az = 0.1 around z = 3? 

The answers to (ii) and (iii) depend on the number of galax- 
ies inside the observed FoV and thus require the apparent 
galaxy positions, unlike in Section |5] where only the appar- 
ent peak flux densities were needed. Therefore this section 
can be regarded as a proto typical a pplication of the S^-SAX- 
Sky simulation (Obreschk ow et al] |2009c). where both radia- 
tive and geometric properties are exploited. 

Throughout the whole section, a "single pointing" refers to 
an observation of a field fixed on the sky with a solid angle 
equal to the instantaneous FoV of the respective observation. 
As in Section |5] such an observation will, in practise, consist 
of several exposures spread over a period longer than the total 
observing time of the pointing. Furthermore we maintain the 
definition that a simulated galaxy is "detected" in a particular 
emission line at n-cr, if the peak flux density of this line 
lies n times above the RMS Jy-noise cr given by eq. ([T]l; re- 
spectively above cr' given by eq. (|8]l in the case of extended 
sources. 

For clarity, the results for the questions (i)-(iii) have been 
collected in a single table (Tab. |5ll, although those questions 
will be explained and discussed separately in Sections 16.11 - 
16.31 The columns of Tab. |5] have again been numbered in 
such a way that this table becomes an extension of Tab. [3] and 
Tab.m i. e. identical columns are given the same column num- 
ber, while new columns are given a new column number. The 
new columns in Tab.|5]are specified as follows. 

(22) Diff'erential number of galaxies at z = 3, detected per 
square degree and unit of redshift at a 3-cr level (or 
higher) in a 24 h single pointing observation. Where no 
value is given (symbol '-'), no galaxy in the S^-SAX- 
Sky simulation can be detected. 
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Emission 


Telescope 


iogl,J ) 




Fraction of MW 


Time to detect 


line 


and band 


[Wm"-] 


[Jy] 


detected in 24 h 


50% of MWs [h] 










3-0- 


10-0- 


3-D- 


10-0- 


(1) 


(2) 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 


HI 


SKAi-LF 


-27.9/ -27.2 


-6.9/ -6.2 


0.0% / 0.0% 


0.0%/ 0.0% 


2.8-10' /I. MO'' 


3.1-10**/ 1-2- 10' 


HI 


SKA2-LF 


-27.9 / - 27.2 


-6.9/ -6.2 


0.0% / 0.0% 


0.0%) / 0.0%o 


Ll-lO*/ 4.4-10'' 


1.2- 10'/ 4.9- 10^ 


CO(l-O) 


SKAi-HF 


-24.0 


-5.2 


5.6% 


0.3% 


360 


4.0-10' 


CO(l-O) 


SKA2-HF 


-24.0 


-5.2 


92.1% 


35.4% 


3.6 


40 


CO(3-2) 


ALMA-3 


-22.6 


-4.3 


48.1% 


7.5% 


25 


280 


CO(4-3) 


ALMA-3 


-22.4 


-4.2 


34.3%) 


5.4% 


43 


480 


CO(5^) 


ALMA-4 


-22.4 


-4.4 


41.5% 


10.2% 


34 


380 


CO(6-5) 


ALMA-5 


-22.8 


-4.8 


16.0% 


4.8% 


360 


4.0- lO'' 


CO(7-6) 


ALMA-5 


-23.4 


-5.4 


9.1%) 


3.5% 


5.3- lO'' 


5.8-10"' 


CO(8-7) 


ALMA-6 


-24.1 


-6.2 


4.6% 


1.1% 


2.8- 10^ 


3.1-10'' 


CO(9-8) 


ALMA-6 


-25.0 


-7.2 


1.7%) 


0.3% 


2.1-10' 


2.3-10** 


CO(10-9) 


ALMA-7 


-26.1 


-8.3 


0.1%) 


0.0% 


8.8-10' 


9.8-10'" 



TABLE 4 

Simulated detectability of the redshifted emission lines emttted by MW progenitors at z = 3. The analysts is based on a complete sample of 173 1 simulated 
MW progenitors. 5 and Sp denote sample medians of the integrated le^ flux and the peak flux density, respectively. Detailed descriptions of each column 

ARE provided AT THE BEGINNE^G OF SeCTIOn[5] FoR HI DETECTIONS TWO VALUES ARE GIVEN, CORRESPO NDING TO THE RAW S'-SAX OUTPUT (/hI = 1) AND THE CORRECTED 

ONE (/hi = 5), AS DESCRIBED IN SeCTIOn IT?) 



(23) Same as ( l22l i but with 10-cr detection limit. 

(24) Absolute number of galaxies in the range z - 2.95-3 .05 
detected at a 3-cr level (or higher) in a 1 yr observa- 
tion using a single pointing for aperture arrays (SKA- 
LF) and 365 distinct 24 h pointings for dishes (ALMA, 
SKA-HF). 

(25) Same as (|24|) but with 10-cr detection limit. 

(26) Number of galaxies with < -22 mag and z = 2.95- 
3.05 inside the instantaneous FoV. If this FoV contains 
less than one object on average, the number is set to 1, 
since we always target known objects. 

(27) Signal-to-noise n of a 24h stacking experiment, where 
the emission lines of all galaxies with Mb < -22 mag 
and z = 2.95-3.05 in the instantaneous FoV are co- 
added. 



6.1. Differential number of galaxy detections in 24 h 

How many galaxies per unit solid angle A and redshift z 
will ALMA and SKA detected in CO and HI emission at 
z = 3 via a 24h-single pointing? Formally speaking, we are 
asking for the differential number count p = dN/(dzdA) at 
z = 3. Due to the discrete number of galaxies, p is computed 
as AA^/(Az AA), where A/V is the integer number of galaxies 
detected within the redshift range Az (around z = 3) and in- 
side the solid angle AA. p is independent of Az and AA if 
three conditions are met: (1) Az is small enough that the dif- 
ferences in the luminosity distances and the effects of cosmic 
evolution can be neglected; (2) AA is large enough to sup- 
press the effects of cosmic variance; (3) the volume spanned 
by Az and AA is large enough that the shot noise on A/V can 
be neglected. We can approximately satisfy these criteria by 
considering the total volume of the S^-SAX-Sky simulation 
contained within the narrow redshift range z = 2.95-3.05. In 
this redshift range the luminosity distance varies by 4% and 
the cosmic look-back time by 82 Mpc within the cosmology 



of the simulation (see Section|2]l. The considered volume ap- 
proximately contains 2.88 ■ 10^ simulated galaxies and covers 
a solid angle of AA ^ 37.2 deg^, which corresponds to the co- 
moving surface area of (500//i Mpc)^ (box size of the Millen- 
nium simulation) that is large enough to suppress the effects 
of cosmic variance. 

p can now be computed for each emission line and tele- 
scope by counting the number of galaxies A/V with peak flux 
densities greater or equal to «-times the Jy-noise of a 24 h- 
observation. For non-resolved sources, this noise cr is cal- 
culated via eq. ^ or, analogously, via cr - (cr VAt)/ V24h, 
where cr VAf is drawn from col. ( fT3] l in Tab. [3] For resolved 
sources, the galaxy-dependent noise cr' of eq. ([8]l needs to be 
adopted instead. The resulting values for p - AN/{Az AA) are 
given col. (l22l i and (|23] | of Tab. |5] Note that these values can 
be computed directly via the SQL-interface of the S^-SAX- 
Sky simulatior0. For example, to get the differential number 
p of CO(3-2) detections using ALMA band 3 (value in Tab.|5l 
col.|22] row 5) the following query can be executed. 



select count(*)/37.2/0.1 
from galaxies Jine 

where zapparent between 2.95 and 3.05 

and coin tfiux_3 *columpeak> 3 *0. 63e-3/sqrt(24 *60) 



Explanation: 37.2 is the FoV AA in deg^ of the S^-SAX-Sky 
simulation at z = 3; O.i is the redshift interval Az; zapparent 
is the apparent redshift of the galaxies including peculiar ve- 
locities; 2.95 and 3.05 are the minimal and maximal values of 
zapparent; cointflux3*columpeak is peak flux density of the 
CO(3-2) line in units of Jy; 3 is the significance level of the 
detection; 0.63e-3 is the value of cryfKt in units of Jy Vmin 
(copied from Tab. [3] col. [13] row 5); 24*60 is the number of 
minutes per day. Note that the values in Tab. |5] may differ by 
up to 30% from those output by the above SQL query, since 
Tab. |5] also accounts for the signal-to-noise decrease in the 
case of extended galaxies (see eq.[8]). 

** http://s-cubed.physics.ox.ac.uk/queries/new?sim=s3_sax 
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(2) 
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(27) 


HI 
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/ 
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HI 
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1 


2 
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1 


21 


CO(3-2) 
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1 


9 
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1 


8 
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ALMA-4 


4.2-10'' 


1.4-10'* 
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58 


1 


12 


CO(6-5) 


ALMA-5 


1.7-10'' 


7.6- 10^ 


50 


22 


1 


6 


CO(7-6) 


ALMA-5 


1.0-10'' 


5.7- 10^ 


22 


12 


1 


3 


CO(8-7) 


ALMA-6 


6.5- 10' 


3.3- 10^ 


11 


5.4 


1 


0.9 


CO(9-8) 


ALMA-6 


4.0- 10^ 


1.7- 10^ 


5.2 


2.3 


1 


0.2 


CO(10-9) 


ALMA-7 


1.5-10^ 


500 


1.6 
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1 


0.03 



TABLE 5 

Simulated detectabujty of the redshifted emission lines of arbitrary galaxies at z = 3. The 1 yr number counts and the stacking analysis depend on the 

RESPECTIVE INSTANTANEOUS FoV. DETAILED DESCRIPTIONS OF EACH COLUMN ARE PROVIDED AT THE BEGINNING OF SeCTI0n|6] FoR HI DETECTIO NS TW O VALUES ARE GIVEN, 
corresponding TO THE RAW S^-SAX OUTPUT (/m = 1) AND THE CORRECTED ONE (/m = 5), AS DESCRIBED IN SeCTIOn IJ^I 



Given the differential number counts p, the absolute num- 
bers of line detections in a 24 h single pointing are obtained 
through multiplying the values of p by the instantaneous 
FoV and by the instantaneous redshift range, which is dictated 
by the instantaneous bandwidth. Hence cols. (l22l i and ( l23T l in 
Tab.|5]cannot be used for a direct comparison of the detection 
rates of the different lines, unless the observed sky field and 
redshift range are smaller than (and hence not limited by) the 
instantaneous FoV and redshift range. The latter case is met, 
for example, when observing a small galaxy group at z - 3. 
In this case, the highest CO detection rates are achieved us- 
ing CO(l-O) [SKA2-HF], followed by CO(3-2) [ALMA-3], 
CO(5^) [ALMA-4], and CO(4-3) [ALMA-3]. 

By comparison, HI detections within a similarly small sky 
field at z = 3 using a 24h SKA observation seem virtually 
impossible. Not a single galaxy of the 2.88 ■ 10^ objects in 
the S^-SAX-Sky simulation at z = 2.95-3.05 has a peak flux 
density above the 3-cr detection limit of SKA2-LF. In fact, 
the 3-cr detection limit corresponds to an HI mass of about 
3 ■ 10'' M0 (assuming an intrinsic line width of 300km s-'), 
which is heavier than the largest HI mass ever ob served in the 
local Universe (e. g. HI Parkes All-Sky Su rvey , Meyer et"al] 
12001 . As we shall demonstrate in Section 16.21 HI nonethe- 
less wins over CO by an appreciable difference during long 
surveys (here 1 yr), where the differences in t he in stantaneous 
FoV become crucial. Furthermore, Section 16.31 reveals that 
even in 24 h observations HI can still be detected at z = 3 
when using SKA in combination with parallel redshift sur- 
veys. 

6.2. Absolute number of galaxy detections in 1 yr 

How many galaxies in a redshift range Az - 0.1 around 
z = 3 will be detected in CO and HI emission during a 1 yr- 
survey? Such a survey can, for example, serve as a mea- 
surement of the angular power spectrum and of the comoving 
space densities fieiC-Z = 3) and D.}\^_(z - 3). 

For such measurements, the precise position of the redshift 
range Az = 0.1 around z = 3 may have to be adjusted in- 
dividually for each emission line to avoid the frequencies of 
radio-frequency interferences (RFIs) and atmospheric absorp- 
tion lines. However, for this analysis we shall use the inter- 
val z = 2.95-3.05 while neglecting RFIs and atmospheric 



absorptions lines. The narrow redshift range Az = 0.1 was 
chosen to isolate this number count analysis from the effects 
of distance variations and cosmic evolution within the sur- 
vey volume. Note, however, that both ALMA and SKA are 
foreseen to yield instantaneous bandwidths corresponding to 
larger ranges in redshift, while easily maintaining our spectral 

resolution of 75 km s-'. 

As revealed in Sections |5] and ISTl SKA-LF is predicted to 
detect virtually no HI at z = 3 within 24 h. On the other hand, 
both ALMA and SKA-HF will be able to detect a significant 
amount of CO. This apparent advantage of dish-based CO de- 
tections over aperture array-based HI detections nonetheless 
vanishes when the instantaneous FoV becomes important. In 
fact, aperture array-based HI searches can perform very long 
exposures of the same, very large sky field, while dishes must 
take many shorter exposures to map a significant sky field in 
the same total observing time. As an example, we here as- 
sume that a 1 yr galaxy survey at z = 3 is performed using 
a 1 yr single pointing for SKA-LF (HI line) and 365 single 
pointings of 24 h each for ALMA and SKA-HF (CO lines). 
In the latter case, the number of detected galaxies is read- 
ily obtained by multiplying the differential number counts 
ANIidzdA) in Tab. |5] (cols. |22l ED by 365, by Az = 0.1, 
and by the instantaneous FoV listed in Tab. [3] (col. fT4l l. In 
the approximation of non-extended sources, the number of 
simulated HI detections in a 1 yr single pointing with SKA- 
LF can be counted by executin g a n ew SQL-query on the S^- 
SAX-Sky database (cf. Section lOT i. For example, for a blind 
HI search using SKA2-LF (Tab. |5] col. 24, row 2) the query 
for this approximation reads 



select count(*)/37. 2*410 
from galaxies Jine 

where zapparent between 2.95 and 3.05 

and hiin tfiux *hilumpeak> 3*0.31 e-3/sqrt(24 *60 *365) 



Explanation: 37.2 is the FoV in deg^ of the S^-SAX-Sky sim- 
ulation at z = 3; 410 is the FoV in deg^ of SKA2-LF; zappar- 
ent is the apparent redshift of the galaxies including peculiar 
velocities; 2.95 and 3.05 are the minimal and maximal values 
of zapparent; hiintllux*hilumpeak is peak flux density of the 
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HI line in units of Jy; 3 is the significance level of the detec- 
tion; 0.31e-3 is the value of cr VAt in units of Jy Vmin (copied 
from Tab. [3] col. [13] row 2); 24*60*365 is the number of min- 
utes per year Note, however, that some values in Tab.|5]differ 
significantly (factor ~ 5) from those output by the above SQL 
query, since Tab. |5] also accounts for the signal-to-noise de- 
crease in the case of extended galaxies (see eq.|8j. 

The results for the absolute number of line detections in the 
range z = 2.95-3.05 are provided in Tab. l5](cols. l24l IBT l. A 
comparison with the differential number counts (cols. l22ll23T l 
highlights the tremendous advantage of aperture arrays. Their 
giant FoV compared to dishes fully compensates the weak- 
ness of HI emission compared to CO emission (e. g. Tab. HI 
col.fTSTl. We further emphasize that the instantaneous FoV of 
SKA-LF is here limited by the computational power of the 
digital back-end. In principle, at least an instantaneous FoV of 
10'*deg^ is conceivable (hemisphere above an elevation of 
30deg), and it seems to be only a matter of time until the 
respective computational resources will become available. 
Therefore, HI surveys with SKA will ultimately be faster than 
any CO survey with ALMA and SKA. 

6.3. Line stacking at z — 3 in 24 hours 

What is the significance of a stacked signal obtained in a 
24 h survey of a redshift range Az = 0.1 around z = 3? Here 
"stacking" refers to the addition of possibly non-detected 
emission lines by using the positions and redshifts of their 
sources drawn from parallel surveys (e. g. optical/infrared). 
Stacking lines of comparable signal and noise ideally in- 
creases the signal-to-noise ratio by V^. It may hence be- 
come possible to measure the summed flux of otherwise non- 
detected emission lines. This technique is illustrated in Fig.|5] 
for three random HI emission lines drawn from the S^-SAX- 
Sky simulation. Stacking typically becomes useful for large 
samples of non-detected lines with redshift uncertainties 6z 
much smaller than the redshift interval spanned by an indi- 
vidual line, i. e. 6z < 10"^ at z = 3 (Fig. |5]). For simplicity, 
we here neglected redshift and position errors, although they 
may be a major problem in real stacking experiments. 

As shown in Fig.|5] a stacked line profile differs from those 
of single emission lines. The maximum flux density of a 
stacked line approximately corresponds to its central flux den- 
sity, since all individual emission lines contribute to its central 
flux, but only the broader lines contribute to its tails. We there- 
fore call a stacked line "detected" if its central flux density is 
detected. Formally, a stacked line composed of individual 
lines is detected at n-cr if 



Perfect line 



Randomly perturbed line 



^/No 



(11) 



where s^ is the central flux density of the i-th line and cr is the 
channel noise (here assumed constant) given by eq. ([U and 
Tab.[3](col.[nii. 

Since stacking experiments require spectroscopic redshift 
measurements from parallel surveys, the latter impose sam- 
ple selection criteria. We here consider a hypothetical stack- 
ing experiment using spectroscopic data from a galaxy survey 
limited to extinction-corrected absolute B-band magnitudes 
Mb < -22. This sample definition approximately matches 
that of the Lyman-break galaxies (LBGs) at z ~ 3, for which 
spectroscopic redshifts were obtained via optica l follow-up 
observations at Keck I and II (ISteidel et al.ll2003h . We stack 
the 24 h emission line signals of all simulated galaxies with 



z = 3.0007 




z = 2.9927 



z = 3.0066 



centered and added 






Fig. 5. — Idea of a line stacking experiment, illustrated with three simulated 
HI emission lines from different galaxies at z a; 3: a low-mass Hl-rich galaxy 
seen edge-on (red), a Hl-rich galaxy seen face-on (green), an intermediate 
mass galaxy seen edge-on (blue). Ideal simulated HI emission lines (left) are 
perturbed with Gaussian channel noise (right). Using the known redshift of 
each individual line (in reality drawn from spectroscopic optical/infrared red- 
shift measurements), the lines are aligned and added to a single line (black). 
The stacked line yields an improved signal-to-noise, roughly a factor V3 
above that of the individual lines. 

Mb < -22 and z - 2.95-3.05, contained in the instantaneous 
FoV of the respective line (see Tab.|3] col.[T3Tl. 

For the CO line observations (ALMA, SKA-HF), the in- 
stantaneous FoV contains, on average, less than one such 
galaxy. In other words, stacking will not allow us to increase 
the signal-to-noise ratio. Only for comparison to HI, we there- 
fore say that the number of stacked galaxies is A^ = 1 and we 
assume that this galaxy has CO fluxes corresponding to the 
geometric average of all simulated galaxies with Mb < -22 
and z = 2.95-3.05. For HI line observations (SKA-LF), the 
instantaneous FoV is very large; i. e. 41 deg^ (SKAi-LF) and 
410 deg^ (SKA2-LF). The requirement of a spectroscopic red- 
shift survey at z = 3 covering such large FoVs, lies far beyond 
current possibilities (e. g. Steidel et al. 2003 considered at to- 
tal sky area of 0.38 deg ). Powerful wide-field spectrographs 
are needed to retrieve spectroscopic data on sky areas com- 
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parable to the instantaneous FoV of the SKA-LF. Even the 
S^-SAX-Sky simulation 'only' covers a FoV of 37.2 deg^ at 
z = 3. In this stacking analysis we therefore linearly extrapo- 
late the number of galaxies and their summed HI fluxes to the 
SKA FoV. 

The number of stacked galaxies and the signal-to-noise 
ratio n of the stacked emission line can be computed using the 
onli ne S QL-interface of the S^-SAX-Sky database (cf. Sec- 
tion |6?T)- For example, for the HI line observed with SKA2- 
LF (Tab. |5] cols. 26, 27, row 2), the respective query reads 



select count(*)/37.2*410as N, 

sum(tabl.hiintHux*tabl.hilumcenter)/(0.31e- 

3/sqrt(24*60))/sqrt(count(*))*sqrt(410/37.2)as n 

from galaxies Jine tab J, galaxies_delu tab2 

where tabl .id-tab2.id 

and tabl.zapparent between 2.95 and 3.05 

and tab2.mag_bdust<-22 



Explanation: 37.2 and 410 respectively denote the FoVs of 
the simulation and the survey; tabl.hiintf}ux*tabl.hilumpeak 
is peak flux density of the HI line in units of Jy; 0.31e-3 is 
the value of cr VAt in units of Jy Vmin (copied from Tab. [3] 
col. [13] row 2); 24*60 is the number of minutes per day; 2.95 
and 3.05 are the minimal and maximal values of the redshift 
tabl.zapparent. This query also calls the table galaxies.delu 
(tab2), which contains the extinction-corrected absolute B- 
band magnitud es mag_bdust calculated in the semi-analyti c 
galaxy model (.Croton et al1l2006l: iDe Lucia & Blaizo3l2007h . 

The significance levels n of the stacked line detections, cal- 
culated via eq. (fTTT l. are listed in Tab.l5](col.l27ll. This analysis 
demonstrates that SKA2-LF can reliably detect HI at z = 3 in a 
24 h pointing using stacking techniques, given large and deep 
spectroscopic redshift surveys. This reveals again the strength 
of the large FoV off'ered by aperture arrays. It should be em- 
phasized that our B-band selection criterion (Mb < -22) ex- 
cludes most of the HI and CO at z = 3. In fact, the stacked 
HI line only traces 2% of the total HI mass in all simulated 
galaxies in the FoV. Larger fractions of the HI mass may be 
studiable by stacking on objects selected by star-formation 
rate indicators, such as used in the wide-area e mission-line 
surveys planned with HETDEX (iHill et al.ll2004 . 

7. DISCUSSION AND CONCLUSION 

7.1. CO detections at z — 3 with ALMA 

Will ALMA meet its primary science goal to detect MW- 
type galaxies at z = 3? Yes it will. Just about. Beginning with 
a semi-analytic model, we selected 1928 simulated galaxies 
that resemble the MW at z = and backtracked their cosmic 
history to a time corresponding to z = 3. In the resulting sam- 
ple of "MW progenitors" or "MW-type galaxies" at z = 3, 
ALMA has roughly a 50% chance of detecting a random ob- 
ject in CO(3-2) emission at a 3-cr level in a 24 h pointing. 
ALMA band 3 is the best choice to achieve this goal, since it 
contains the CO(3-2) line and since it is the only band con- 
taining another low order CO transition at z = 3, i. e. CO(4- 
3). If the instantaneous bandwidth is split into two windows, 
covering CO(3-2) and CO(4-3), and if these two lines are co- 
added ("single source stacking", see Section |5]l, the odds of 
detecting a MW in 24 h can be increased to 60%. These pre- 
dictions remain similar if instead of using a model for MW- 
progenitors at z = 3, we simply imagined the actual MW at a 



cosmological distance equivalent to z = 3, because the total 
H2 mass only changes by about 20% between the two cases 
(see Tab. [T]). Whether those predictions will be met signifi- 
cantly depends on the final, currently uncertain sensitivity of 
the ALMA receivers and on the actual transparency of the at- 
mosphere at the different bands. 

As expected from its small instantaneous FoV, ALMA is a 
relatively slow survey instrument. On average, only about 1-2 
general galaxies (not just MWs) per day will be detected in a 
blind CO survey between z = 2.95 and z = 3.05 (from col. 24 
in Tab. |5]l. To make effective use of ALMA as a CO survey 
instrument, it is hence crucial to preselect a sample using CO 
indicators, such as tracers of star formation. 

7.2. C0( 1-0) and HI detections at z — 3 with SKA 

SKA-HF - should its frequency domain be extended up to 
28.8 GHz - will provide a unique way of detecting CO at z = 
3. In fact SKA-HF searches for CO(l-O) promise to become 
much more effective in terms of number of detected objects 
than ALMA searches for CO(3-2). This result is based on 
the assumption that only the core of SKA-HF is used to keep 
most CO(l-O) sources non-resolved, and it already accounts 
for the fact that SKA-HF yields a low antenna efficiency 
(see col. 12, Tab.O at 28.8 GHz due to its limited dish surface 
accuracy. The power of SKA-HF compared to ALMA relies 
in two main reasons: firstly, SKA-HF will ultimately reach 
a total collecting area three orders of magnitude larger than 
that of ALMA, and secondly the instantaneous FoV (i.e. the 
beam size) is an order of magnitude larger for CO(l-O) than 
for CO(3-2). 

On the other hand, HI detections at z = 3 using SKA-LF are 
a mixed blessing. On the downside. Section |5] revealed that 
SKA-LF will be virtually unable to pick up a MW progenitor 
at z - 3. Even when using the full SKA2-LF array and when 
assuming that the typical MW progenitor yields 5-times more 
HI than predicted by the S^-SAX model, it would still take 5 
years (4.4 ■ 10"^ h, see col. 20, Tab. |4|i of effective observing 
time to pick up such a MW progenitor at only 3-sigma signif- 
icance. We can conclude that MW progenitor studies at z = 3 
in HI will remain virtually impossible using the SKA, while 
CO detections with ALMA seem possible within 24 h. On the 
other hand. Section |6] suggests that blind searches for HI in 
general galaxies (not just MW progenitors) at z = 3 will be 
more effective than any conceivable CO survey with ALMA 
and SKA-HF. In fact, the huge instantaneous FoV of SKA- 
LF, which is roughly 6 orders of magnitude larger than that of 
ALMA, compensates for the low fluxes of HI lines compared 
to those of CO lines. As a result SKA2-LF promises to detect 
above 10^, perhaps even above 10^, galaxies atz = 2.95-3.05 
in a 1 year blind search for HI, while ALMA will find less than 
10^ CO emitters in the same time. Furthermore, if the posi- 
tions and redshifts of the most Hl-rich galaxies in the SKA-LF 
FoV are already known from a parallel spectroscopic survey, 
then stacking techniques can be applied to detect HI in less 
than 24 h (col. 27,Tab.|5]l. 

Finally, it is worth considering the implications of extend- 
ing the mid-frequency array SKA-MF down to 355 MHz to 
observe the HI line at z = 3. In this case, the wave front of 
the HI line can be sampled completely, implying a geome- 
try factor close to = 1 compared to eg =0.1 of SKA-LF 
(see Tab. O. The implied gain in sensitivity is roughly com- 
pensated by the 10-times smaller collecting area of SKA-MF 
compared to SKA2-LF. However, SKA-MF would still have 
better point source sensitivity because of its lower receiver 



14 



temperatures (Tree ~ 50 K, see Tab.|2]i and because it will pro- 
vide a better image quality due to reduced side-lobes above 
the horizon. The instantaneous FoV is also likely to be signif- 
icantly larger due to the smaller station size. In conclusion, 
SKA-MF promises interesting features for observing HI at 

7.3. Closing words 

This joint study of ALMA and SKA has highlighted re- 
markable differences in their performances and applicabili- 
ties. In general, ALMA is most powerful at targeted ob- 
servations, while SKA will be very suitable for long blind 
searches, thus placing ALMA preferentially in the surround- 
ings of galaxy evolution, while SKA will reach deep into large 
scale cosmology. Thus ALMA and SKA will be heavily syn- 
ergetic. Together those instruments will set a supreme long- 
term standard in speed and sensitivity over the whole atmo- 
spherically transparent frequency domain between 70 MHz 



and 920 GHz; and together they will resolve many outstand- 
ing questions across all redshifts of the star-forming universe. 

With minor exceptions, this paper has been restricted to 
z = 3 and to the current benchmarks for ALMA and SKA. 
However, the methods and on-line simulation tools presented 
here can be applied to any other redshift and any other tele- 
scope. This paper therefore sets the stage for the statistical 
comparison of future cold gas surveys with simulations. Such 
a comparison is a crucial, if not the only, way to verify physi- 
cal theories of galaxy formation against the empirical reality. 
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